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We experimentally study the effects of the anisotropic Rydberg-interaction on D-state Rydberg 
polaritons slowly propagating through a cold atomic sample. We observe the interaction-induced 
dephasing of Rydberg polaritons at very low photon input rates into the medium. We develop a 
model combining the propagation of the two-photon wavefunction through our system with nonper- 
turbative calculations of the anisotropic Rydberg-interaction to show that the observed effect can 
be attributed to pairwise interaction of individual Rydberg polaritons at distances larger than the 
Rydberg blockade. 


Long-range and spatially anisotropic dipole-dipole 
(DD) interactions enable new approaches for prepar¬ 
ing and exploring strongly correlated quantum systems 
[1]. Magnetic DD interaction couples individual nuclear 
spins to nitrogen-vacancy centers in diamond [2, 3] and 
is observed in dipolar gases of ultracold atoms [4, 5]. 
Electric DD interaction determines the long-range in¬ 
teraction between polar molecules [6] or Rydberg atoms 
[7, 8] and may allow to investigate phenomena such as 
quantum magnetism [9-11] and topological phases [12] 
in these systems. Recently the angular dependence of 
the DD interaction between single Rydberg atoms has 
been fully mapped [13]. Here we study, for the first 
time, the anisotropic DD interaction between slowly 
travelling polaritons coupled to a Rydberg D-state via 
electromagnetically-induced transparency (EIT) [14]. 

Rydberg-EIT has emerged as a powerful approach to 
realizing few-photon optical nonlinearities [15-18]. This 
novel technique enables a variety of optical quantum 
information applications such as highly efficient single¬ 
photon generation [19], entanglement generation between 
light and atomic excitations [20] , single-photon all-optical 
switches [21] and transistors [22, 23], and interaction- 
induced photon phase shifts [24]. Additionally, it al¬ 
lows us to probe novel phenomena such as attractive 
interaction between single photons [25], crystallization 
of photons [26], or photonic scattering resonances [27]. 
The electric DD interaction between pairs of Rydberg 
atoms has been studied extensively in the perturba¬ 
tive van-der-Waals regime [28-30]. Rydberg-EIT exper¬ 
iments have so far mostly employed Rydberg S-states, 
where the interaction has only very weak angular de¬ 
pendence. Recent experiments simultaneously prepare 
Rydberg atoms in S- and P -states [31] or two different 
S'-states [22, 23, 32]. The interaction between these en¬ 
ergetically well-separated states enables novel entangle¬ 
ment schemes for atomic systems [33] and increased flex¬ 
ibility in the manipulation of weak light fields [34, 35]. 

In this work, we employ Rydberg EIT to investigate 
the anisotropic interaction between individual Rydberg 
polaritons. We observe that D-state Rydberg polaritons 
are decoupled from the EIT control field due to their 
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FIG. 1. (a) EIT level scheme with week probe field £ p 

and strong control field D. (b) Frequency scan of the weak 
probe field S p showing the EIT transmission window aver¬ 
aged over the full pulse time. The transmission shows a 
strong nonlinearity with incident photon rate, (c) Geomet¬ 
ric scheme of our setup consisting of EIT lasers at 780 nm 
(wq « 6 fi m) and 480 nm (wq ~ 12 fin i) focused onto a thermal 
cloud (30 (i K) of 87 Rb atoms with extensions of a z = 80 fi m 
and (j r = 25 fim. The grey shadow shows the anisotropic 
blockade region caused by a |100D 5 / 2 ) Rydberg excitation, 
(d) Time dependent transmission on two-photon resonance 
for different probe photon input. 


pairwise interaction. We show that this effect is signifi¬ 
cant for incident probe photon rates corresponding to a 
mean spacing between polaritons larger than the Ryd¬ 
berg blockade distance. We extract the scaling of the 
decoupling rate with EIT control Rabi frequency and 
principal quantum number of the Rydberg D-state. Fi¬ 
nally, we present a model which treats the complicated 
pair interaction as a dephasing process. Combining de¬ 
phasing rates extracted from nonperturbative pair poten¬ 
tial calculations with numerical propagation of the two- 
polariton wavefunction we model our experimental obser¬ 
vations. Our EIT scheme consists of a few photon probe 
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field £ p and a strong coupling field Cl resonantly coupling 
the levels illustrated in Fig. 1(a). Resulting transmission 
spectra of £ p for a frequency scan of the probe laser over 
EIT resonance employing the I 8 OD 5 / 2 ) Rydberg state for 
two different photon rates R{ n are shown in Fig. 1(b). For 
Ri n = 0.15photons/yus we observe high transmission and 
a narrow linewidth, from which we extract a decoherence 
rate 7 ^/ 2 ^ = 200 kHz, originating from the movement of 
the atoms due to the finite temperature and stray electric 
and magnetic fields over the cloud. For higher probe in¬ 
put rates the EIT transmission decreases, due to the self¬ 
blockade of propagating polaritons as observed in exper¬ 
iments using Rydberg 5-states [16, 18, 19]. In contrast, 
here we observe a qualitatively new effect present only 
when we apply the EIT scheme to Rydberg D- states, 
namely a time dependence of the transmission on EIT 
resonance (Fig. 1(d)). For very low photon rates on the 
order of Ri n = 0.15 photons//is, where we rarely have 
two photons inside the medium, we measure a constant 
transmission, which is determined by Cl and the decoher¬ 
ence rate j gr . If we increase R { n , we observe a decay of 
transmission over time which gets faster with increasing 
rate Ri n . We attribute this effect to interaction-induced 
coupling to degenerate Zeeman sublevels which leads to 
polaritons being converted to stationary Rydberg exci¬ 
tations inside the cloud. These impurities shift the Ry¬ 
dberg levels of the surrounding atoms, preventing other 
polaritons from propagating through the cloud [21-23]. 

To investigate the dependence of the decay of trans¬ 
mission on the photon rate Ri n , we calculate an effective 
optical depth ( OD e ff ) of our medium on EIT resonance 
(A = 0) by taking the logarithm of the transmission 
(Fig. 2(a)). To account for the three effects described 
above, we write OD e ff as 

OD e ff = ODdec + OD n i(R in ) + ODdph(Rin,t ), (1) 

where the terms ODdec , OD n i and ODdph represent the 
contributions of decoherence, blockade-induced nonlin¬ 
earity, and dephasing. In the following we investigate 
the last contribution. Neglecting saturation effects for 
the highest values of it^ n , we approximate the increase 
in OD due to dephasing to be linear in time and write 
ODdph = Rod T, with Rod as creation rate of optical 
density by decoupled impurities. 

In Fig. 2(b) the extracted rates Rod are plotted ver¬ 
sus Ri n for measurements with different control field 
Rabi frequencies Cl coupling to the \88D 5 / 2 ,mj = 5/2) 
state. For measurements with CI/2tt = 16.6 MHz and 
CI/2tt = 26.3 MHz, respectively, Rod scales quadrati- 
cally with R{ n over the probed range. For Cl/2 tt = 
6.1MHz and Cl/2ir = 10.8 MHz we see deviations from 
the quadratic dependence for Ri n exceeding 1.5 and 2.7, 
respectively. Considering the delay time in the medium 
given by T de i ay = °7 7e (where 7 e /27T = 6.1MHz is the 
decay rate of intermediate state I 5 P 3 / 2 )) and finite initial 
EIT transmission at t = 0, the onset of these deviations 
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FIG. 2. (a) Effective optical density of the \88D 5 / 2 ) EIT 

feature for different photon rates at fixed CI/2tt = 8.3 MHz. 
The solid lines are linear fits to the data to extract the rates 
Rod of additional OP e // due to the creation of impurities, 
(b) Dependence of Ro d on Ri n for different Rabi frequencies 
Cl. Parabolic fits (solid lines) determine the rate constant 
C(Q) for each dataset. The vertical dashed lines correspond to 
mean photon numbers inside the medium exceeding 2 where 
the quadratic dependence breaks down, (c) Rate constants C 
for different Cl. For the different principal quantum numbers 
n we observe the same scaling according to Eq. (3) with k = 
1.67(4). Dashed lines are fits to the data. 


seems to coincide with rates corresponding to a mean 
number of photons in the medium exceeding 2 (vertical 
lines in Fig. 2 (b)). This quadratic dependence of Rod cm 
Ri n suggests that the observed dephasing is a two-body 
effect. In particular, we will later relate the dephasing 
to the probability I'lpddl 2 of finding two polaritons in the 
Rydberg state at the same time for distances larger then 
the blockade radius 77 , given by 

I 1 M 2 = • (2) 

Because of this dependence, we introduce a rate constant 
C(Cl), relating Rod and Ri n in the quadratic regime via 
Rod — C(Cl) • Rf n , which we obtain from fits to our data 
(Fig. 2(b)). The Cl dependence of the extracted rate con¬ 
stants is shown in Fig. 2(c) for measurements with prin¬ 
cipal quantum numbers n = 80, n = 88 and n = 100 . 
To compare the results for different principal quantum 
numbers, we extract the dependence on Cl by a fit of the 
form 


c(n) = a-n~ k , (3) 

which yields k = 1.67(4) for all the different n. How- 
ever, for the prefactor a we observe a strong scaling with 
n, indicating significantly larger dephasing for larger n. 
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In order to connect our observations to the anisotropy 
of the Rydberg interaction we calculate Rydberg pair 
potentials through diagonalization of the electrostatic 
DD interaction Hamiltonian [36]. To reduce computa¬ 
tion time we diagonalize the Hamiltonian in the coordi¬ 
nate system aligned with the atomic separation, where 
the total magnetic moment M = mnj\ + mj 2 is con¬ 
served. Calculated pair potentials for |8RD 5 / 2 ; 8 CLD 5 / 2 ) 
and |IOOD 5 / 2 ; 1001 ) 5 / 2 ) are shown in Fig. 3(a) and (b). 
The anisotropic DD interaction couples states with dif¬ 
ferent magnetic quantum numbers in the D-pair state 
manifold, resulting in the observed splitting of the man¬ 
ifold in addition to the overall level shift. When con¬ 
sidering the coupling of a laser field with given direction 
and polarization one has to account for the electric dipole 
selection rules by calculating the overlap of the coupled 
pair state in the fixed coordinate system defined by the 
direction of light propagation, with the eigenstates in the 
presence of interaction in the interatomic coordinate sys¬ 
tem [29]. The colored shadow of the potentials shows 
the overlap of the \nD^/ 2 ^mj = 5 / 2 ;nT 5 / 2 ,Tnj = 5/2) 
pair state in the fixed coordinate system, rotated into the 
interatomic frame via Wigner d-matrices, with the new 
eigenstates for an angle 6 = 60° between the interatomic 
axis and the light propagation direction. It can be seen 
that the coupled state is projected onto multiple new 
eigenstates. Unlike for isotropic 5-states this projection 
depends strongly on the angle 0 [29]. One consequence 
of the anisotropic interaction is that the blockade volume 
for D-states is not spherical as the interaction potential 
can no longer be described by a single Cq coefficient. Sec¬ 
ondly, the coupled pair state is not an eigenstate in the 
presence of interaction and will thus experience a time 
evolution under the influence of the interaction Hamilto¬ 
nian. 

To incorporate the interaction induced dynamics with 
polariton propagation under EIT, we reduce the full in¬ 
teraction potential to an effective level structure as shown 
in Fig. 4(a). In the weak-probe limit, it is sufficient 
to consider only two photons simultaneously inside the 
medium [18, 27]. Accounting for level shifts and selection 
rules, we determine the anisotropic blockade distance 
rb(0) where the laser-coupled pair state is shifted out 
of the EIT bandwidth (Fig. 4(b)). From this we obtain 
an effective interaction potential between two slow light 
polaritons V(z,r±). This term results in the blockade- 
induced nonlinearity of the medium. Outside the block¬ 
aded region, where both the overall shift and the split¬ 
ting of the D 5/2 Zeeman-manifold are smaller than the 
EIT bandwidth, we describe the evolution of the cou¬ 
pled pair-state by an effective dephasing rate T(z,r±) 
of the Rydberg polaritons into localized Rydberg exci¬ 
tations. Here, z = z\ — Z 2 and r_ j_ = r_ l 5 i — r_ i _ j2 are 
the relative coordinate along and perpendicular to the 
propagation direction. To obtain T(z,r±) we calculate 
the time evolution of a stationary Rydberg pair, in the 
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FIG. 3. Calculated pair state potentials of (a) 

|80D 5 / 2 ; 8 OD 5 / 2 ) and (b) |100D 5 / 2 ; 100D 5/2 ) pair states (gray 
lines). The color shading represents the projection of the 
|raji ,2 = 5/2, 5/2) state with respect to the the light propa¬ 
gation direction onto the new eigenstates in the presence of 
interaction for an angle of 0 — 60° between the interatomic 
axis and the coupling laser beam. 


initial state coupled by the fixed control field, at given 
distance and angle between light propagation direction 
and interatomic axis from the full interaction potential. 
Although this dynamics is fully coherent, the revival of 
the initial population appears only on time scales long 
compared to the polariton propagation time due to the 
large number of states in the D-state manifold. Hence, 
the initial time evolution on experimentally relevant time 
scales is well described by a spatially varying dephasing 
rate. The important result is that for D-states T(z,r_\_) 
is large at distances well beyond the blockade volume 
(Fig. 4(b)). In contrast, the same approach results in 
vanishingly small dephasing rates for 5-states. 

For the propagation dynamics we solve numerically the 
full set of propagation equations for the two-body wave 
function [18, 25], including the dephasing rate T(z,rj_) 
as a decay of the amplitude ^dd{ r 1 , 1 * 2 ) of two Rydberg 
excitations. We assume a homogeneous distribution of 
atoms inside the finite-size medium of length L = 4 cr z 
and include imperfect single-photon transmission due to 
the decoherence of the 2 -photon transition. Neglect¬ 
ing probe-beam diffraction due to the interaction, we 
solve polariton dynamics (Fig. 4(c)) for different r± with 
Vid(z) = V(z, r±) followed by averaging over the r± dis¬ 
tribution determined by the Gaussian transverse profile 
with waist w e jj = 7 jim (corresponding to the waist aver¬ 
aged over length L) of the probe beam. Then, the rate 
of events J\f that at least one photon is converted into 
stationary Rydberg excitation is proportional to the am- 
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FIG. 4. (a) We model the anisotropic DD interaction be¬ 

tween two polaritons as an effective potential V(r = {z,r±}) 
capturing the interaction-induced blockade effect, and a posi¬ 
tion dependent dephasing rate T (r) out of the pair-state cou¬ 
pled by the EIT control laser, (b) Example plot showing the 
anisotropic blockaded volume and the dephasing rate outside 
this region for n = 80 and Q/2tt = 25 MHz. (c) Numerical 
simulations showing: For z > 0 probability distribution asso¬ 
ciated with two Rydberg excitations \^dd(z, R)\ 2 . For z < 0 
the product T(z)\ipdd(z 1 R)| 2 , where R — (zi + z^)/2 is the 
center of mass coordinate. Both quantities are presented in 
arbitrary units for r± = 4.2 fim, ft/2jr = 12 MHz, n — 80 and 
probe-beam waist w e ff— 7 [im. (d) Comparison of C(Q) from 
numerical simulations (circles) with fits to experimental data 
from Fig. 2(c) (dashed lines). 


plitude of the two-polariton wavefunction i/jdd- 

AT = f dridr 2 T(ri - r 2 )|V> dd (ri,r 2 )| 2 . (4) 

Jv 

We normalize the Rydberg wave function ipdd with the 
incoming photon flux as in Eq. (2). A single dephasing 
event increases on average the optical depth within the 
medium by OD irn for the incoming photons. Therefore, 
the change in transmission by a single dephasing event is 
exp (~OD dec - OD sa t(R in )) (l - e~ ODim ). This reduc- 
tion appears with the rate A/", and therefore the initial 
time evolution for the averaged transmission behaves as 

T(t) = e -° D dec-OD sat {R in ) exp (i _ e -OD irn ' jj 

and leads to a rate constant C( ft) = 

AT (l — e~ ODirn ) /R? n . Furthermore, we observe, 

that a finite life time of the Rydberg impurity results in 
an effective saturation of the transmission. However, the 
full time evolution for the transmission including higher 
number of excited Rydberg impurities is extremely chal¬ 
lenging and beyond the scope of the present manuscript. 
These calculations provide important insights into 
the behavior of the dephasing. First, the averaged 
optical thickness OZA m > 1 of a dephasing event is 


sufficient to strongly block the medium. Therefore, the 
decay mainly follows the probability to absorb at least 
one impurity. Second, there appears a characteristic 
distance with the highest probability for the excitation 
of an impurity Rydberg state given by the competition 
between the higher dephasing rates at shorter distances 
and the suppression to find two Rydberg excitations 
due to the blockade effect, Fig. 4(b-c). The latter is 
strongly affected by the polariton dynamics inside the 
medium, as has been previously discussed in terms of 
a diffusive behavior [18]: at the entrance of the two 
photons into the medium, the probability to find two 
Rydberg excitations is purely determined by the block¬ 
ade due to interactions 'ipdd(z) ~ 1/(1 — xYid(z)) where 
hx = — i( 7 e /fl 2 + l/ 7 e ) [27]. This dip in probability 
broadens during propagation due to correction to the 
linear behavior of the polariton dispersion, while the 
single polariton losses provide an overall reduction of 
the amplitude (see Fig. 4(c)). These effects strongly 
depend on ft. In addition, photons inside the medium 
are compressed due to the slow light velocity, which con¬ 
tributes an additional factor fl -4 (Eq. (2)) to the scaling 
of C(ft) with ft. Both described effects combined explain 
the numerical results presented in Fig. 4(d). We find 
qualitative agreement between theory and experiment 
without any fitting parameters. While for low n and 
small values of ft the agreement is excellent, for larger ft 
we observe a discrepancy, and moreover, the theory does 
not reproduce the strong scaling with main principal 
quantum number n measured in the experiment. We 
expect the reason for this discrepancy to be the fact that 
for large ft and n the AC-stark shift and broadening 
of the Rydberg lines, caused by coupling to the 5 P 3/2 
manifold by the control field, becomes comparable to the 
splitting between the nD 5 / 2 and nD 3 / 2 manifold which 
scales with n~ 3 . In this case, our two-step approach 
of first calculating the interaction potentials and then 
incorporating them into the polariton propagation does 
not capture the full evolution of the system. 

In summary, we have investigated the effect of 
anisotropic DD interaction between Rydberg D-state po¬ 
laritons. Interaction-induced decoupling into stationary 
Rydberg excitations results in a decrease of transmission 
on EIT resonance over time. This effect is relevant to 
all Rydberg experiments employing non-P-states [31] or 
Forster resonances [23, 32], where the anisotropy of the 
Rydberg interaction will always result in coupling to ad¬ 
ditional levels. Our theoretical approach to include full 
Rydberg pair state potentials in numerical two-photon 
propagation in the form of an effective potential and 
anisotropic decay rate yields qualitative agreement with 
our measurements and thus is a useful tool for treat¬ 
ing the complicated interaction between Rydberg polari¬ 
ton and stationary excitations in general. The fact that 
we observe the interaction-induced state-mixing on the 
few-photon level is a promising result for experiments 
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on Rydberg-dressing [10, 11] and engineered polariton- 
interaction [37] using Rydberg states with orbital angu¬ 
lar momentum. More detailed study of the anisotropic 
coupling will be possible in storage and retrieval exper¬ 
iments [21, 38], enabling control over number and posi¬ 
tion of stored excitations. In this scenario, it becomes 
particularly interesting to employ echo techniques [39] to 
probe the coherent spin evolution of interacting Rydberg- 
polaritons. 
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